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Abstract The Nonlinear stability of triangular equi- 
librium points has been discussed in the generalised 
photogravitational restricted three body problem with 
Poynting-Robertson drag. The problem is generalised 
in the sense that smaller primary is supposed to be an 
oblate spheroid. The bigger primary is considered as 
radiating. We have performed first and second order 
normalization of the Hamiltonian of the problem. We 
have applied KAM theorem to examine the condition of 
non-linear stability. We have found three critical mass 
ratios. Finally we conclude that triangular points are 
stable in the nonlinear sense except three critical mass 
ratios at which KAM theorem fails. 

Keywords Nonlinear Stability: Triangular Points: 
Generalised Photogravitational:RTBP:P-R Drag 



1 Introduction 

The simplest form of the three-body problem is called 
the restricted three-body problem(RTBP), in which a 
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particle of infinitesimal mass moves in the gravitational 
field of two massive bodies orbiting according to the 
exact solution of the two-body problem. In the circu- 
lar problem, the two finite masses are fixed in a coor- 
dinate system rotating at the orbital angular velocity, 
with the origin (axis of rotation) at the centre of mass 
of the two bodies. Lagrange showed that in this ro- 
tating frame there are five stationary points at which 
the massless particle would remain fixed if placed there. 
There are three such points lying on the line connecting 
the two finite masses: one between the masses and one 
outside each of the masses. The other two stationary 
points, called the triangular points, are located equidis- 
tant from the two finite masses at a distance equal to 
the finite mass separation they are stable in classical 
case. The two masses and the triangular stationary 
points are thus located at the vertices of equilateral tri- 
angles in the plane of the circular orbit. There is a 
group of enthusiasts who want to setup a colony at L5 
point of the Earth-Moon system. As already noted, be- 
cause L4 and L5 are the stable points of equilibrium, 
they have been proposed for sites of large self-contained 
"Space colonies" , an id ea developed and advocated by 
the late lO 'Neil] (|l974h . The three body problem have 
an interesting application for artificial satellites and fu- 
ture space colonization. Triangular points of the Sun- 
Jupiter or Sun-Earth system would be convenient sites 
to locate future space colonies. Application of results 
to realistic actual problem is obvious. 

The classical restricted three body problem is gen- 
eralized to include the force of radiation pressure, the 
Poynting-Robertson(P-R) effect and oblateness effect. 
The photogravitational restricted three body prob- 
lem arises from the classical problem when at least 
one of the interacting bodies exerts radiation pres- 
sure, for example, binary star systems(both primaries 
radiating). The photogravitational restricted three 
body problem under different aspects was studied by 



ladzievskii ( 1950 ). Chernikov J 1970h . ?. Schuermanl of La in the restricted three body problem . Hence we 



1980), llshwar and Kushvahl (2006). iKushvah Sharma and I ahwatej study nonlinear stability of triangular points in 



2007al) 



The Poynting-Robertson drag named after John 
Henry Poynting and Howard Percy Robertson, is a pro- 
cess by which solar radiation causes dust grains in a 
solar system to slowly spiral inward. Poynting ( 19031 ) 
considered the effect of the absorption and subsequent 
re-emission of sunlight by small isolated part icles in the 
solar s ystem. His work was later modified bv lRobertson 
(1937|) who used precise relativistic treatments of the 
first order in the ratio of the velocity of the particle to 
that of light. 

The location and stability of the five Lagrangian 
equilibrium points in the planar, cir cular res t ricted 
three-body problem was investigated by Murray ( 19941 ) 
when the third body is acted on by a variety of drag 
forces. The approximate locations of the displaced equi- 
librium points are calculated for small mass ratios and 
a simple criterion for their linear stability is derived. 
They showed if a\ and 03 denote the coefficients of the 
linear and cubic terms in the characteristic equation 
derived from a linear stability analysis, then an equilib- 
rium point is asymptotically stable provided < a\ < 
CI3. In cases where a\ is approximately equal to or a± 
is approximately equal to 03 the point is unstable but 
there is a difference in the e-folding time scales of the 
shifted L4 and L5 points such that the L4 point, if it 
exists, is less unstable than the L5 point. The results 
are applied to a number of general and specific drag 
forces. They have shown that, contrary to intuition, 
certain drag forces produce asymptotic stability of the 
displaced triangular equili brium points, L4 and L5. 
Ishwar and Kushva h (2006) examined the linear stabil- 
ity of triangular equilibrium points in the generalised 
photogravitational restricted three body problem with 
Poynting-Robertson drag and conclude that the trian- 
gular equilibrium points are unstable due to Poynting - 
Robertson drag . Kushvah Sharma and Ishwarl (|2007bl ) 
performed higher order normalizations in the general- 
ized photogravitational restricted three body problem 
wi th Poynting-Robertson drag. 

Deprit and Deprit-Bartholomei ( 1967 ) investigated 



the nonlinear stability of triangular points by applying 
Moser's modified vers i on of Arnold's thcorcm(1961). 
Bhatnagar and Hallanl (|l983h studied the effect of per- 
turbati ons on the nonlinear stability of triangular 
points. Ilshwarl (| 19971 ) studied nonlinear stability in the 
generalized restricted three body problem. His problem 
is generalized in the sense that the infinitesimal body 
and one of the primaries have been taken as o blate 
spheroid. Subba Rao and Krishan Sharmal (1997) ex- 
amined effect of oblateness on the non-linear stability 



our problem. 

To Examine the nonlinear stability of triangu- 
lar points we use d the KAM theoremjthe work of 
Kolmogorov ( 1957 ) extended by Arnold ( 1961 ) , Moser ^ 



(1962)]. Vloser's conditions are utilised in this study by 
employing the iterative scheme of Henrard for trans- 
forming the Hamiltonian to the Birkhoff's normal form 
with the help of double D' Alembert's series. We have 
found the second order coefficients in the frequencies. 
For this we have obtained the partial differential equa- 
tions which are satisfied by the third order homoge- 
neous components of the fourth order part of Hamil- 
tonian H4 and second order polynomials in the fre- 
quencies. We have found the coefficients of sine and 
cosine in the homogeneous components of order three. 
They are critical terms. We have eliminated these crit- 
ical terms by choosing properly the coefficients in the 
polynomials. Then we have obtained the values of the 
coefficients A, B, C occurring in the fourth order part 
of the normalized Hamiltonian in KAM theorem. Wc 
have applied KAM theorem to examine the conditions 
of nonlinear stability. Using the first condition of the 
theorem, we have found two critical mass ratios /i c i , /i C 2 
where this condition fails. By taking the second order 
coefficients, we have calculated the determinant D oc- 
curring in the second condition of the theorem. From 
this, we have found the third critical mass ratio /i C 3 
where the second condition of the theorem fails. We 
conclude that triangular points are stable for all mass 
ratios in the range of stability except three critical mass 
ratios where KAM theorem fails. The stability condi- 
tions are different from classical case and others, due 
to radiation pressure, oblateness and P-R drag. 



2 First Order Normalization 



We used IWhittakerl (|1965f ) method for the transforma- 
tion of H2 into the normal form 



Equations of motion are as in llshwar and Kushvah 
(2006) and given by 
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+ fi) +y 2 ,r 2 2 = (x + fi-lY+y 2 , 
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(x + fJ,)[(x + n)x + ; 



+ x — ny, 



y[(x + u)x + yy] 



+ y + n(x + a) 



W = (l-M)(l-gi) = ma_ 

masses of the primaries, A 2 = 



< |, mi,77i2 be the 
t 2 - be the oblate- 
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ness coefficient, r e and r p be the equatorial and polar 
radii respectively r be the distance between primaries, 
Cd = 299792458 be the dimcnsionlcss velocity of light, 

/ F \ 

qi = (1 — -pf j be the mass reduction factor expressed 
in terms of the particle's radius a, density p and radia- 
tion pressure efficiency factor \ (hi the C.G.S. system) 
i.e., q\ = 1 — 5,6 — -. Assumption q\ = constant 
is equivalent to neglecting fluctuation in the beam of 
solar radiation, the effect of the planet's shadow, obvi- 
ously q± < 1. Triangular equilibrium points are given 
by U x = 0, U y = 0, y ^ 0, then we have 



x < 1 



Wi[(l-n)(l + %A a ) +m(1 



3^(1 - ^)y £o 




and the Hamiltonian is H = —L-\-p x x-\-p v y, where 
PxiPy are the momenta coordinates given by 



Px 



Py 
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^ 1( 5 2 [2 M - 1 - M (l - Mi)f + 7(1 - M )#] 
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For simplicity we suppose q\ = 1 — e, with |e| << 1 
then coordinates of triangular equilibrium point L4 can 
be written in the form 
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that, we change x — » a;* + a; and y 
a = + /i, 6 = y* so that 



.</* 




A 2 + 



2A 2 e 



(9 + 7 V-^ l} 
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Expanding L in power series of x and y, we get 



(5) 



L = L Q + L x + L 2 + L 3 + ■ ■ 
H = H + H 1 +H 2 + H 3 + 
= -L+p x x + p y y 



(7) 
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where 7=1 — 2/i. We shift the origin to L4. For 

y. Let 
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wher e = 4, — yn = ±8(1 — ir) ^ an< 4 ^ = <7i ^ 3 , 
as in " 
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The Lagrangian function of the problem can be writ- 
ten as 

1 n 2 
L = -^{x 2 + y 2 ) + n(xy- xy) + — (x 2 + y 2 ) 

! (l~M)gi + ii + ^ 2 
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where L , Li,L 2 ,L 3 
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43e 
36 
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(14) 



(i 2 + y 2 ) , . . \ n 2 , 2 2 , 
2~ hn(xy-x?/j + y(a; + y) 

-Ex 2 - Fy 2 - Gxy (15) 



L 3 = -gj {x 3 T! + 3x 2 yT 2 + 3xy 2 T 3 + y 3 T 4 + 6T 5 } 

(16) 

£4 = -^{^i^ + ^^y + GiVaa; 2 ?/ 2 

+4^ 4 ^ 3 + 24^ 6 } (17) 



where 



E 



16 



2 - 6e - 3A 2 - 



31A 2 e _ mW x 

2 ~~ 6V3 



+7{2e + 1M2 + ^ + ^,] (18) 



G 



-1 
16 

7 {6e 
V3 



10-2e + 2L4 2 - 



717 A 2 e 67Wi 



18 
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The second order part H 2 of the corresponding 
Hamiltonian takes the form 



Hi 



Px t ^ +<yPx -xp y )+ Ex 2 + Fy 2 + Gxy (21) 



To investigate the stability of the motion, as in 



Whittakerl (|1965T ) , we consider the following set of linear 



equations in the variables x, y: 



9H 2 
dx ' 



-X Py = 


dH 2 
dy ' 


i.e. 






AX = 
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X 




X = 


y 

Px 

. Py . 





X.r 



dH 2 
dH 2 



(22) 



2E G A -n 

G 2F n A 

-A n 1 

—n -A 1 
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Clearly \ A\ = 0, implies that the characteristic equa- 
tion corresponding to Hamiltonian H 2 is given by 

\ i + 2(E+F+n 2 )\ 2 +AEF-G 2 +n A -2n 2 {E+F) = 

(24) 

This is characteristic equation whose discriminant is 

D = A(E + F + n 2 ) 2 -4{4:EF-G 2 + n 4 -2n 2 (E + F)} 

(25) 

Stability is assured only when D > 0. i.c 

(i < Mc D - 0.221896e + 2. 103887^2 + 
0.493433eA 2 + 0.704139 + 0.401 154eWi (26) 

where fi Co = 0. 03852 l,(Routh's critical mass ratio) 
When D > the roots ±iuj 1 and ±iuj 2 (cJi, w 2 be- 
ing the long/short -periodic frequencies) are related to 
each other as 



2 1 2 
+ u) 2 



je 3 7 A 2 83eA 2 Wi 
2 + 2 + 12 24\/3 



(27) 



27 27 7 2 9e 9 7 e 117 7 A 2 



+ 



16 16 8 8 16 

Tj, Nj, (i = 1, . . . , 5, j = 1, . . . , 6) are as in lKushvah Sharma and Ishwaij 241e ^ 2 + 3514/1 _ 55 ^ W i 
( 2007bh . 32 16^3 

(0 < w 2 < -p < wi < 1) 
v2 



16 



(28) 
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From and ( 
satisfy 



it may be noted that Uj (J = 1, 2) 



2 _ 4e _ 107eA 2 2 7 e 2514^ 
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27 

16 32e 
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729 + 243^3/ 
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81 



1880e^ 2 _ 2720W/i \ 
729 ~ 2187V3/ 



(29) 



Alternatively, it can also be seen that if u = wiw 2 , then 
( [281) gives 



7 



1 + 



4e 107eA 2 



9 



25^i 



+ 



27 27V3 
1579eA 2 557^1 
324 qTT 

32e , 208A 2 



27 243 



81 



1880eA 2 320WA 2 
729 + 243^37 



(30) 



Following the method for re ducing H 2 to the normal 
form, as in IWhittakerl ljl965h ,use the transformation 



X = JT 



where X 



(31) 



x 

y 

Px 

Py 



,J=[J, 



y'Jl<ij<4! 



T 



Qi 

Q2 
Pi 

P-2 



J21 \ / sin 0i + J 22 < / sin 2 

V ^1 v ^2 

+^23^2/1^1 cos 0i + J24V2I2UI2 sin0a(35) 



3 Second Order Normalization 



In order to perform Birkhoff's normalization, we 



where Jy are as in lKushvah Sharma and Ishwarl(|2007bh .| 

P, = (2/^) 1/2 cos0 z , Q t = (^) 1 /2 s i n ^ ) (i = 1,2) 

The transformation changes the second order part of 
the Hamiltonian into the normal form 



H2 — Ullli — U2I2 



(32) 



The general solution of the corresponding equations of 
motion are 



const. 



±Wj + const. 



1,2) 



(33) 



If the oscillations about L4 are exactly linear, the 
Eq. ( I33|) represent the integrals of motion and the cor- 
responding orbits will be given by 



Jl 3 \/2wi/i COS 0i + Juy/2u)2h COS0 2 (34) 



use H enrard's method (jDeprit and Deprit-Bartholome 
( 1967t )) for which the coordinates (x, y) of infinitesi- 
mal body, to be expanded in double D'Alembert series 



Er 



B 



1,0 



y 



J2n>i Bn where the homoge- 



neous components B^' and B®' 1 of degree n are of the 
form 

^2 I 1^ _I 2 ^2 [ C n-m,m,p,q COS (p0i + # 2 ) 
0<m<n (p><?) 

sin (p<Pi + qfa)] (36) 

The conditions in double summation are (i) p runs over 
those integers in the interval < p < n — m that 
have the same parity as n — m (ii) q runs over those 
integers in the interval — m < q < m that have the 
same parity as m. Here I\, I 2 are the action momenta 
coordinates which are to be taken as constants of in- 
teger, 0i, 2 are angle coordinates to be determined 
as linear functions of time in such a way that 0i = 

w l + Y,n>\ f2n(h,h), 02 = -0>2 +Y,n>l 92n(h,h) 

where u>i, uj 2 are the basic frequencies, f 2n and g% n are 
of the form 



f%n — ^ f 2(n-m),2m^l "^2™ (37) 
0<m<n 

92n = ^ 9 2{n-m),2m,^\ ™^2™ (38) 



0<m<n 

The first order components B\'° and B®' 1 are the val- 
ues of x and y given by ( [M]) ( [35]). In order to find 
out the second order components B 2 ' , B 2 ' 1 we consider 
Lagrange's equations of motion 

d ,dL. dL „ d ,dL. dL _ , . 

^-^ = ' 5^)-^ = ° (39) 



i.e. 



S — 2ny + (2P — n 2 )ir. + 
x + 2ni + (2F — n 2 )y + G.t 



9^3 i 9Lj 
dx dx 



9« cty 



(40) 



Since x and y are double D'Alembert series, x 3 x (J > 
0, /c > 0, j + k > 0) and the time derivatives x,y, x,y 



6 



are also double D'Alembert series. We can write 



where x, y, x, y are homogeneous components of degree 

. r l/2 r l/2 . 

n m L-y , ±2 i.e. 



E 



^M^l + / 2 + / 4 + 



rj>l L 



-W 2 + 92 + 54 + 



(41) 



We write three components Xi,X2, £3 of x 



Xi 



X-2 



•f'3 



Ull- 



dBi 



1,0 



w 2 - 



dBi 



1.0 



CJl- 



UJ 2 - 



DB. 



1,0 



dBl>° 

-— w 2 



U(p 2 

8Bl'° dB{>° 

r /2-^7 — 



32" 



DB^ + h 8 -^- 

d(j>i 



1,0 



52- 



1,0 



902 



where 



D = uj 1 



8 



(jj 2 - 



a 



(901 U(p 2 

Similarly three components xx,x 2 , £3 of x are 
xi = D 2 B{'°, x 2 = D 2 Bl'°, 



(44) 



(45) 



n2 D i,o , o f 9 B 1 d B 1 

x 3 = D B 3 ' + 2^1/2^73 2u; 2 ff2 » , 2 

002 



50? 



In similar manner we can write the components of y, y. 
Putting the values of cc, y, x, y, x and y in terms of dou- 
ble D'Alembert series in Ea. (HP]) we get 



D* + 2E - 1 - -A 2 ) B L 2 



2\1+1m) D-C4\BV 



X 2 (46) 



3 

r 



2 ( 1 + -A- 2 ] : ' Gj //_!" 

3 
2' 



+ ( D 2 + 2F - 1 - ^A 2 ) B2' 1 = Y 2 (47) 



X 



where 
dL 3 



dx 



Yo 



- B i ,v—b 1 



8L 3 



dy 



D l,0 D 0.1 

x=B t ' ,y=B^ 



These are two simultaneous partial differential equa- 
tions in B^' and B 2 ' ■ We solve these equations to 
find the values of B^° and B° A , from (@BJ> and (@7|> 



AiA 2 i?2' = $2, A1A2-B 



0,1 

2 



-\&9 



(48) 



D 2 + a;. 



where Ai=D 2 +w 2 ,A 2 
$2 = (-D 2 + 2F- n 2 )X 2 + (2nD - G)Y 2 



VPs 



(2nD + G)X 2 - (D 2 + 2E- n 2 )Y 2 



(49) 
(50) 



The Eq.( l4"5|) can be solved for B\'° and B®' 1 by putting 
the formula 



1 



(42) 


AiA 2 


(43) 


where 


1,0 


A p ,g = 




1 



A 



p.q 



cos(p0i + q<j> 2 ) 
or 

sin(p0i + g0 2 ) 



[w 2 - (wip - uj 2 q) 2 ] [oj 2 - (viP ~ ^2<?) 2 ] 

provided A p , q 7^ 0. Since A^o = 0, A ,i = 
the terms cos 0i , sin 0i , cos 2 , sin cf>2 are the critical 
terms. $2 and ^2 are free from such terms. By condi- 
tion^) of Moser's theorem kiuii + k 2 0J2 7^ for all 
pairs (ki,k 2 ) of integers such that |fci| + \k 2 \ < 4, 
therefore each of Wi,W2,Ui ± 2o7 2 ,w 2 ± 2uji is differ- 
ent from zero and consequently none of the divisors 
Ao.o, Ao,2j A2,o ; A11, Ai _i is zero. The second order 
components B^^B®' 1 are as follows: 



Bo 



B. 



0.1 



Tih + r 2 h + r$h cos 20i 
+r 4 / 2 cos 202 + r b ll ,2 ll /2 cos(0i - 2 ) 
+r & l{ /2 ll /2 cos(0i + 02 ) + r 7 h sin 20i 
+r s I 2 sin 20 2 + rg/f 2 /^ sin(0i - 2 ) 
+ri O /i 1/2 / 2 1/2 sin(0i + 2 ) (51) 

-|si/i + s 2 I 2 + S3A cos 201 

+ S 4 ^2 COS 202 + S5/i /2 / 2 1/2 COS(01 - 2 ) 

+s 6 /i /2 /2 1/2 cos(0i + 02 ) + s 7 h sin 20i 
+s a h sin20 2 + sgl^ 2 /^ 2 sin(0i - 2 ) 
+ S i O /i 1/2 /2 1/2 sin(0i + 2 )} (52) 



where r ,; , s$, (i = 1, . . . , 10) are as in Kushvah Sharma and Ishw arj 
Using transformation a; = Bl' + B2' and 
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y = B®' 1 + B®' 1 the third order part H 3 of the Hamil- 
tonian in l\^ 2 , l^ 2 is of the form 



H 3 = A 3fi lV 2 +A 2tl I 1 I^+A lj2 I^h+A , 3 I? z (53) 

We can verify that in Ea.f IS*3l A 3 n vanishes inde - 
pendently as in Deprit and Deprit-Bartholomel ( 1967 ). 
Similarly the other coefficients A 2A , Ai j2 , j4o,3 are also 
found to be zero independently. 



r l/2 



1/2, 



7-3/2 



4 Second Order Coefficients in the Frequencies 

In order to find out the second order coefficients 
f 2> o, /o,2> 9 2 Qj <7n 2 m the polynomials f? and w e have 
done as in lDeprit and Deprit-Bartholomel (|l967l ). Pro- 
ceeding as ( HHJ , we find 



A 1 A 2J B 3 1 '° = $ 3 -2/ 2 P-2 ff2 g 

A 1 A 2 B°' 1 = 9 3 -2f 2 U-2g 2 V 
where 

$ 3 = \d 2 + 2F- n 2 ] X 3 + \(2nD - G)l Y 3 



2{nD + G) 



X 3 + 



L> 2 + 2nE - n 2 



(54) 



(55) 



(56) 



Y 3 (57) 



[D 2 + 2F- n 2 ] 
+ {2nD - G) 



2ni.o 



LUl- 



d 2 Bj_ 



dB { ; 



0,1 



d 2 B"> L dB L { 
uji „ h n- 



d0? 



(58) 



[D 2 + 2F - ra 2 ] 



2 R 1,0 



W 2 



crB\_ 



dB° { 
1 dfa 



(2nD - G) 



u 2 - 



d 2 B^ 



dBi 



1.0 



d 2 B] fi 8B" 
u>i — Trm n 



(59) 



U = -(2nD + G) 
+ [D 2 + 2E-n 2 ] 



o.i 



d<j>\ 

8 2 B°' 1 dB\> 
9^1 901 



(60) 



r 



(2riL> + G) 



,2 D 1,0 



^'2 



c# 2 



0^ 



0.1 



- [L> 2 + 2£ - n 2 ] 



2 7,0,1 



uj 2 - 



d 2 B]_ 



dB{ 



1.0 



d(j) 2 



(61) 



dy 



(62) 



i.e. 



*3 



2 2 



y x + T 2 a;& + y 2/ + -g-a; 
iV 3 2 iV 4 3 , dT 5 dN 6 

+T y + 

T 2 , T 4 2 ^2 3 ^3 2 

+T 3 a;2/ + y2/ + ~q x + V 



(63) 



iV 4 



iVi 



+-?ran/ z + — y A + — h 
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dy dy 



(64) 



(|S"6"|) and (fST)) are the partial differential equations 
which are satisfied by the third order components 



-Bg'^-Bg' 1 and the second order polynomials / 2 ,g 2 in 
the frequencies. We do not require to find out the 
components B\' a and S3' 1 . We find the coefficients 
of COS <pi , sin (j>i , COS (j> 2 and sin (j> 2 in the right hand 
sides of ( [56l. ([5T)l . They are the critical terms , since 
A10 = Ao,i = 0. We eliminate these terms by choos- 
ing properly the coefficients in the polynomials 



/2 = /2,0-fl + fo,2h, 92 = 52, O^l + 30,2^2 

Further, we find that 

(coefficient of COS0i in <& 3 ) 



fo,; 



/2,0 



52,o 



50,2 



2( coefficient of COS0i in P) 
(coefficient of COS0 2 in $3) 
2( coefficient of COS0 2 in Q) 
(coefficient of COS0 2 in ^3) 
2( coefficient of cos0 2 in Q) 



(65) 

A (66) 
B (67) 
G (68) 



where 



B 



C 



where A\ t i,B\ t i and Gi^, (i 
pendix I 



5 Stability 



A 1A + (Ai >2 + A x $y)e + (A 1A 
+A h5l )A 2 + {A h6 + Ai j7 7)Wi (69) 

Bi A + {B h2 + B h3l )e+(B 1A 
+B h5l )A 2 + (Bi i6 + 51,77)^1 (70) 

G 1 , 1 + (G 1)2 + Gi,37)e+(C M 
+Gi, 57 )A 2 + (Gi, 6 + Ci,ry)Wi (71) 

7) are as in Ap- 



1. 



The condition(i) of KAM theorem fails when u>i 
and io\ = 3w 2 



2uj 2 
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5.1 Case(i) 

When uj\ = 2tt>2 
Then from ( l?2"| and ( 



we have 



„ , _27 3e 117^2 _ 221WA 
5 " 4 ~ 2 ~ 4 ~ 15V3 J 



+ 



4 100 100 
I 407e 12A 2 2399111', 



5.2 Case(ii) 

When oji = 3u> 2 
Proceeding as f !5.1[) . we have 

M ^ 4 2 4 - 20 J 

(27 93e 29794a 119\/3~wV 
+A< ^ 4 ~ 100 100 + 10 

9 393e _ 27A 2 _ 4777^1 
"100 200 ~ 100 ~ 400%/3 

Solving for [i-, we have 



(72) 



727 107e 3021 A 2 4291 Wi\ 



120^5^ 

(73) 



25 200 25 200V3 
Solving for /i we have 

Mci = 0.024294 -0.312692e 
-0.036851A 2 + 1.001052Wi (74) 



(75) 



/i c2 = 0.013516 - 0.29724e 
-0.019383A 2 + 1.007682Wi (77) 

Normalized Hamiltonian up to fourth order is 

H = wiJi - ^ 2 / 2 + -(AI 2 + 2Bhh + C7f ) + . . . (78) 

Calculating the determinant D occurring in condition 
(ii) of KAM theorem, we have 

D = -(Awl + 25wicj 2 + CujI) 

Putting the values of A, B and C and if u = u)\uj 2 , we 
have 

644m 4 - 541u 2 + 36 , «. 
J= 8(4, 2 ~l)(25, 2 -4) + ^ + ^ 
+(^4 + D 5l )A 2 + (D 6 + D r j)W! (79) 

The second condition of KAM theorem is satisfied if , 
in the interval < /i < /i c o, [where ^co as in (|26p ] 



the mass parameter does not take the value fic3, which 
makes D = 0. To find /i C 3, we note that when e = 
A2 = W\ = 0, then from ( l?9"]) . D becomes zero if and 
only if 

644u 4 - 541u 2 + 36 = 
This implies that 



541 - V199945 
1288 



Writing u 2 = 27 ^ 1 16 7 ^ , 7 = 1 — 2/i and then solving 
above, we get 

7 = 7o = 0.978173 n = fio = 0.010914 . . . (81) 

When e,A 2 ,Wi are not zero, we assume that 13 is zero 
if 



0.072863 = u (say) (80) 



fi = fi + a\e + a 2 ^4 2 + a 3 Wi 

7 = 70 - 2(aie + a 2 ^4 2 + aaWi) 
where 70 = 1 — 2^0 

u 2 = u + («i + aiM 2 )e 



(82) 



(83) 



(76) with 



+ (u 3 + a 2 u 4 )A 2 + (u s + a z UQ)Wi (84) 



27 2 , 9 ,9 _ 117(1 -7o 2 ; 



Ul = 16 7 ° + 8 7 ° + °' U3 



16 



u 2 = u 4 



277o 



27to 
4V3' 



u 5 



27 7o 2 + 165 7o + 35 
16V3 



and cti, (i = 1, 2, 3) are to be determined. From ( PT9"|) . 
Z? is zero when 

644u 4 - 541u 2 + 36 , „ N 
= 8(4 U 2 -1)(25 U 2 - 4) +(Z?2 + 2337)6 
+(D 4 + Z? 57 )A 2 + (£> 6 + £77)^1 = (85) 

Making use of ( in ( l8"S"|) and equating to zero the 
coefficients of e, A 2 and Wi, we get 



1 



-{(1288u -541)ui 



m 2 (1288u -541) 
+8{D° 2 + £>°7 )(4tt - 1)(25« - 4)} (86) 
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ft: 2 



0(3 



- , 4 (1288l-541) {^ 1288 -°- 541 ^ 3 
+8{D A + D°7o)(4 Uo - l)(25«o - 4)} (87) 



^{(1288^-541)^ 

m 6 (1z88u — 541) I 

+8{D° 6 + D°7o)(4ix - l)(25«o - 4)} 



where (n = 2, 3, 4, . . . , 7) are D„ given as in Ap- 
pendix II, as evaluated for the unperturbed problem. 
Numerical computation yields, 

ai = -0.120489 . . . , a 2 = -0.373118 . . . , 

a 3 = 2.904291 . . . 

Then we have 

Mc3 = Mo + ftie + "2^2 + a 3 Wi 

= 0.010914 -0.120489e 

-0.3731 18^2 + 2.90429 lWi (89) 

Hence in the interval < /i < /i c o, both the condi- 
tions of KAM theorem are satisfied and therefore the 
triangular point is stable except for three mass ratios 
[i ci {i = 1,2,3). 



6 Analytical Study 

6.1 Observation I 

Consider A 2 = 0,qi = 1, (W\ = 0) then problem re- 
duced to the classical restricted three body problem. 
From equation ( rj| ( [5]) we get 

1 ^V3 
x = 2 - M, V = ±~2 

from ( I26p stability is assured when fi < /i c q where 
(J-cO = 0.038521. The relation between u>x,cj2 in 
(I27) ) (l25 )l are given by 



u\+ul = l, ^1 = ^1-^ (01)) 



(0 < uj 2 < —/= < UJl < 1) 
v2 

From (EH), (ED) (155)) we nave found that the tri- 
angular points are stable in the range of linear stability 
except the three mass ratios 



A*ci 

Mc2 
A*c3 



0.024294 
0.013516 
0.010914 



(91) 
(92) 
(93) 



and the D occurring in the second condition of KAM 
theorem we have found from ( I79p 



D 



644). 



54k 



36 



8(4m 2 -1)(25m 2 -4) 



(94) 



where u = uj\uj 2 

All the a bove results, are exactly similar with the 
results as in lDeprit and Deprit-Bartholomd (|1967l ). 

Now we have e = 1 — q\,W\ = ^7~, suppose 
D = Dq. We draw the figure ( [1]) which describes the 
instability range in classical case and figure ( [5]) views 
the points u> x = 0.924270, uj 2 = 0.381742, D Q = , 
when A 2 = 0, qi = 1 and i < ui < 1, < uo 2 < 
the value of 7 = 0.978173. 

6.2 Observation II 

Consider the case when A 2 = 0, q\ 7^ 1, (Wi 7^ 0) i.e. 
photogravitational restricted three body problem with 
P-R drag when bigger primary is supposed to be radiat- 
ing body and small primary is being spherical symmet- 
ric. The coordinates of triangular equilibrium points 
are given by 



x 



1 



wM(i- M ) + /4] 

3/x(l - n)x y Q 



(95) 



WxS 1 ^- 1 



y = yo V- Hi-M 

this result coincides with ISchuerman ( 198dl ). where 



(96) 



x 



- Mj y 
m, y = 



Substituting e 



1 



±5(1 

1 - = 

p = Mci, (» = 0, 1, 2, 3),A 2 = in ( El, ( E3 
( I89[) . we have found that the triangular equilibrium 
points are stable in the range of stability except three 
mass ratios 



Mc2 



^c3 



0.024294- 0.312692(1 - q x ) 
0.976732(1 -q x ) 



0.013516- 0.29724(1 - qi ) 
0.994062(1 - qi ) 



+ 



0.010914- 0.120489(1 - qi ] 
2.87259(1 - qi) 



+ 



(97) 



(98) 



(99) 



We have observed from table U) and figure ([3]), the 
mass ratio increases, accordingly as the radiation pres- 
sure increa ses, these results a re similar but not identical 
to those of Papadakisl ( 1999 ). 
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. 6 



Fig. 1.— A 2 = 0, qi = 1 





Fig. 3. — Stability region /i, c .;, {i = 1, 2, 3) — <?i , when 
A 2 = 0, W x i= 

6.3 Observation III 

When A 2 7^ 0, <ji = 1, (Wi = 0) i.e. in this obser- 
vation we have considered the smaller primary as an 
oblate spheroid, the radiation pressure(P-R drag) is not 
considered. The triangular equilibrium points are given 

by 



1 - 2/i - A 2 



± 



V3 
2 



(100) 
(101) 



which are similar but not identical to results as in 



Bhat nagar and Hallanl (|1983f ) and lChandra and Kumar 
(|2004f ) . In this case triangular equilibrium points are 
stable in the nonlinear sense except three mass ratios 
at which Moser's condition fails. Which are given by 



Mc3 



0.024294-0.036851A 2 (102) 
0.013516-0.019383A 2 (103) 
0.010914-0.373118^2 (104) 



The stability region are shown in the diagram A 2 — 
Hci{i = 1,2,3), ( [3]), the outer line is correspond- 
ing to /j, c i, second line due to fi c2 and innermost 
line is due to /x C 3 it is clear from table (0) the /x 

decreases as A 2 increases. These results agree with 

Markellos Papadakis and Perdios (1996); Bhatnagar and Hallanfr 
( 19831) 



Fig. 2.— A 2 = 0,gi = l,ui = 0.924270, uo 2 = 
0.381742, D = 



6.4 Observation IV 

When A 2 ^ 0, qi ^ 1{W\ ^ 0) this is the most gen- 
eralized case which is being considered. The triangular 
equilibrium points are given by ( 2|) , ( clearly they 
are the functions of oblateness coefficient A 2 and P-R 
drag term Wi . 
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0.2 0.3 0.4 0.5 



Fig. 4. — Stability region /U c j, (i = 1,2,3) — A 2 , when 
q x = = 



in the nonlinear sense in the range of linear stabil- 
ity except for three mass ratios // c ,, (i = 1,2,3) at 
which KAM the orem fai ls , whi le they are unstable in 
linear case [see iMurravl (|1994l ); llshwar and Kushvah 
( 2006hl. These results agree with those found by 



Deprit and Deprit-Bartholomei ( 1967 ) and others 
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Substituting e = 1 — qi,Wi 



(i-M)e 



H = fj, ci , (i 



0, 1, 2, 3) in ([731), ([771) (EH), we get the new formulae 



Mc2 



Mc3 



0.024294- 0.312692(1 - qi ) 

0.976732(1 - qi ) 



- 0.036851A, 



0.013516- 0.29724(1 - qi ) 

0.994062(1 - qi ) 



- 0.019383A 2 + 



Cd 



(105) 
(106) 



0.010914- 0.120489(1 - qi ) 
0.373118A 2+ 2 - 87259 ( 1 -^) (107) 

Cd 



Using ( HD5|l-( [TOTJ we have drawn /z - A 2 - qt, 3D 
diagrams ([5]). You can see in the first diagram, the 
uppermost plane is due to /i c i, middle plane is due to 
\i c i and innermost plane is due to /x C 3 , second view value 
of /ico = .035829. From these diagrams, we reached at 
the conclusion that the stability region is reduced due 
to P-R drag and oblateness effect of smaller primary. 
But still the triangular equilibrium points are stable in 
the range of linear stability except three mass ratios 
at which KAM t heorem fails, w h ile they are unstable 
in line ar case [see IMurravl (|1994T ) ; llshwar and Kushvah 
(120061) ]. 



7 Conclusion 



Using IWhittakerl (j 19651 ) method we have seen that the 
second order part H2 of the Hamiltonian is transformed 

into the normal form H2 = — LO2J2 and the third 

1 12 1/2 

order part H3 of the Hamiltonian in l x ,I 2 ' zero. 
We conclude that the stability region is reduced due 
to P-R drag and oblateness effect of smaller primary 
But still the triangular equilibrium points are stable 
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Table 1 A 2 = 0, Ql £ 1, (Wi ^ 0) 


qi 


Mel 


Me2 


Mc3 


0.95 


0.00866 


-0.001346 


0.00488921 


0.96 


0.011786 


0.0016263 


0.006094 


0.97 


0.014913 


0.0045987 


0.007299 


0.98 


0.018040 


0.0075712 


0.008504 


0.99 


0.02117 


0.010544 


0.00970878 


1.00 


0.024294 


0.013516 


0.0109137 



Table 2 A 2 ^ 0, gi = 1, (Wi = 0) 



A 2 


Mel 


Mc2 


Mc3 


0.0 


0.024294 


0.01352 


0.010914 


0.1 


0.020609 


0.01158 


-0.026398 


0.2 


0.016924 


0.009639 


-0.06371 


0.3 


0.013239 


0.007701 


-0.101022 


0.4 


0.009554 


0.005763 


-0.138334 


0.5 


0.005869 


0.003825 


-0.175645 


0.6 


0.002184 


0.001886 


-0.212957 


0.7 


-0.001501 


-0.000052 


-0.250269 




Fig. 5. — Both the graphs show the stability region fx c i, (i = 1, 2, 3) — q\ — A 2 , second graph view fi c o = .035829 
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Appendix I 

Coefficients Ai :i} Bi yi and C\^, (i = 1 . . . , 7) are given by: 

9 259w 2 

A l.l = - , ; , n 9W -, , r 9N + 



8(-l + 2^ 2 ) 2 (-l + 5^ 2 ) 24(-l + 2w 2 ) 2 (-l + 5w 2 ) 
205w? 31w? 



18(-l + 2w 2 ) 2 (-l + 5w 2 ) 18(-l + 2w 2 ) 2 (-l + 5a; 2 ) 
1 13w 2 



A 1 



2 



36(1 - 2co 2 ) 2 (-l + 5tul) 18(1 - 2co 2 ) 2 (-l + 5w 2 ) 



13wf 167wf 



+ 



27(l-2w 2 ) 2 (~l + 5o> 2 ) 72(l-2co 2 ) 2 (-l + 5u; 2 ) 

W7uf 

108(1 - 2w 2 ) 2 (-l + 5w 2 ) 

421w 2 



2(1 - 2^ 2 )3(-l + 5^ 2 ) 2 32(1 - 2co 2 )3(-l + 5w 2 ) 2 
19wf 8141559w? 



2(1 - 2cj 2 ) 3 (-1 + 5^ 2 ) 2 32(1 - 2w 2 ) 3 (-l + 5w 2 ) 2 
29wf 



^1,4 = 



(l-2w 2 ) 3 (-l + 5w 2 ) 2 

1319 12639w 2 



436(1 - 2w 2 ) 2 (-l + 5ujf) 436(1 - 2w 2 ) 2 (-l + 5w 2 ) 
14275cjf 799cj6 



436(1 - 2co 2 ) 2 (-l + bul) 218(1 - 2w 2 ) 2 (-l + 5tof) 
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57 525i^ 2 

1,5 52(1 -2w 2 ) 3 (-l + 5w 2 ) 2 52(l-2co 2 )3(-l + 5u; 2 ) 2 
475w? 1559w? 



26(1-2cj 2 ) 3 (-1 + 5cj 2 ) 2 26(1 - 2cj 2 ) 3 (-1 + 5w 2 ) 2 
283w? 



Ai.e 



13(l-2w 2 )3(-l + 5w 2 ) 2 

2747^ 2 , 41(9 + 4luI) 



10368V3(-1 + 2cj 2 ) 9216V3(-l + 2w 2 ) 2 
93899(9 + 4w 2 ) 12875w 2 (9 + 4w 2 ) 



331776V3(-l + 2^ 2 ) 82944V3(-1 + 2w 2 )- 



1337 779wi(9 + 4^ 2 ) 41(9 + 4w 2 ) 

- - e1 A57 i i n « + moco AT/' i i n ox, + ~ " ~ 



6144\/3(-l + 2w 2 ) 10368V3(-l + 2w 2 ) 2 18432V3(-1 + 2co 2 ) 2 

227347o> 2 (9 + 37259(9 + 4a; 2 ) 

331776V3(-l + 2w 2 ) 82944\/3(-l + 2w 2 ) 
6517w 2 (9 + 4w 2 ) 



3072V3(-1 + 2cj 2 ) 2 (4w 2 - w 2 ) 

43ct>iW2 



6(1 - 5w 2 )(-l + 2w 2 )(l - 5cj 2 )(1 - 2w 2 ) 

32w 3 u;3 



-61,2 



3((1 - 5w 2 ))(-l + 2w 2 )(l - 5w 2 )(l - 2w 2 ) 
309wio;2 



8(1 - 5co 2 )(l - 2w 2 )(l - 5w 2 )(l - 2w 2 ) 

5904wiW2 407w 3 w 3 . 



Bl,3 



(-1 + 2w 2 )(9 + 4w 2 ) 2 6(1 - 5tul)(l - 2tul)(l - 5w 2 )(l - 2w 2 ) 
1800^2 



Bl,5 



(-l + 2^ 2 )(9 + 4^ 2 ) 2 

10083 - 61407(k; 2 u; 2 + 400800a^| - 3035216^f - 260802u;?a;f 
+ 8^x^(9 - 59w 2 + 62cof + 40^6) (9 - 59w 2 + 62w| + 40wf ) 

247^1^2 

1,4 "~ 4(1 - 5w 2 )(l - 2a; 2 )(l - 5w 2 )(l - 2^ 2 ) 

6817^^3 

36(1 - 5^ 2 )(1 - 2^ 2 )(1 - 5w 2 )(l - 2w 2 ) 

1800wiw 2 



(-l + 2^ 2 )(9 + 4^ 2 ) 2 

, -89211 + 2042998w 2 w| + 1028577^|w 2 + 16052098w?w| + 1215804a;fa;| 



32wiW2(-l + 5^ 2 ) 2 (-l + 5w 2 ) 2 (-9 + 14w 2 + 8wJ)(-9 - 14w 2 + 8w|) 
1599V3(9 + 192cj 1 w 2 +w 2 ) 



Bifi 



512cj 1 w 2 (-1 + 2cj 2 )(9 + 4w 2 ) 2 
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3^3(2398599 - 90316804 - 369i^4 + 574^| + 15744o^| + 3284 
1,7 ~ 51244(-l + 24)(-9 + 144 + 84) 

192(-41601 + 414)4 



512w?w£(-l + 2 w 2)(-9 + 144 + 84) 



~ 9 205o; 2 2 



Cl,3 



Cl,7 



8(-l + 2w|) 2 (-l + 5a;|) 24(-l + 2w|) 2 (-l + 5w|) 
2054 + 314 



18(-l + 24) 2 (-l + 54) 18(-l + 24) 2 (-l + 54; 



134 



36(1 - 24) 2 (-l + 54) 18(1 - 24) 2 (-l + 5^ 
134 1674 

4 



27(1 - 24) 2 (-l + 54) 72(1 - 24) 2 (-l + 54) 

1074 

108(1 -24) 2 (-l + 54) 



4214 



2(1 - 24)3(-l + 54) 2 32(1 - 24)3(-l + 54)- 
194 4074 



2(1 - 24) 3 (-l + 54)2 16(1 - 24) 3 (-l + 54)2 



294 

-f 



(l-24)3(-l + 54) 2 

1319 126394 



436(1 - 24) 2 (-l + 54) 436(1 - 24) 2 (-l + 54) 
142754 7994 



436(1 - 24) 2 (-l + 54) 218(1 - 24) 2 (-l + 54) 



Cl,5 — =777- _ „-„, - ; _ 9N „ + 



52(1 -24) 3 (-l + 54)2 52(1 -24) 3 (-l + 54)2 
4754 15594 



26(1- 24) 3 (-l + 54)2 26(1- 24) 3 (-l + 54)2 

2834 

+ 13(1-24)3(^1 + 54)2 



287V3(-3 + 324 + 484 ) 
1,6 ' 10244(-9 + 144 + 84 



V382^ 2 (3 _ 384 + 164 + 964) 

512(9 + 44) (-4 + 44) (4 - 24)2 

, 3V34(-142911 + 1951104 + 747284 + 667844) 



512(9 + 44) (-4 + 44) (4 - 24)2 
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Appendix II 



567(-151 + 16^ 2 ) f 1620864 2507364 

2 ' 16384(-1 + 2^)2(9 + 4^2)2 + 884736 { (-1 + 2^ 2 ) 2 (1 - 2u\) 
706482 71663616000 



(-l + 2w 2 ) 2 (4u; 2 -w 2 ) (-l + 2c^)(9 + 4o;2) 
8062156800 1074954240 



+ 
+ 



(-l + 2^ 2 )(9 + 4^ 2 ) 2 (-l + 2w 2 ) 2 (9 + 4w 2 ) 

112969617408 

(1 - 5^ 2 ) 2 (-l + 2w 2 ) 2 (9 + 4w 2 )(l - 5w|) 2 (-l + 2^)2(9 + 4^2) 

17146183680 
(9 - 14w 2 - 8w: 



1028577wf 



16(1 - 5co 2 ) 2 (-l + 2w 2 ) 2 (9 + 4luD(1 - 5cj 2 ) 2 (-1 + 2w 2 ) 2 (9 + 4w 2 ) 

8026049^ 

"8(1 - 5w 2 ) 2 (-l + 2^ 2 ) 2 (9 + 4^ 2 )(1 - 5^ 2 ) 2 (-l + 2^ 2 ) 2 (9 + 4w 2 ) 



303951^? 



, .8 
U> 2 



4(1 - 5w 2 ) 2 (-l + 2^ 2 ) 2 (9 + 4c^ 2 )(l - 5^ 2 ) 2 (-l + 2^ 2 ) 2 (9 + 4w 2 ) 

D 3 r sow 

3 8192(-1 + 2w 2 ) I (2 Wl +w 2 )(u;i + 2w 2 )(9 + 4^ 2 ) 2 

6883328 i 

~ (1 - 5w?)(-l + 2w 2 )(9 + 4^ 2 )(1 - 5w 2 )(-l + 2u; 2 )(9 + 4o> 2 ) J 

f 706240 289737 . 2 _ 2 . 530841600 
+ 147456 I (-1 + 2^ 2 ) + (l-2c^ 2 ) ( 1 ^ (-1 + 2 W 2 )(9 + 4w 2 ) 2 
59719680 59719680a; 2 . 3317760 



(-l + 2^ 2 )(9 + 4w 2 ) 2 (-l + 2w 2 )(9 + 4w 2 ) 2 (-l + 2o> 2 ) 2 (9 + 4w 2 ) 
| 71516160 24772608 



9-14w 2 -8w| (w 2 -4^ 2 )(-l + 2a; 2 )(9 + 4^ 2 ) 
22637076480 



'(1- 


5 W 2 )(- 


-1 + 2cj 2 )(9 + 4w 2 )(1-5w 2 )(- 
100200^w| 


-l + 2w 2 )(9 + 4w 2 )/ 


(1- 


5c 2 )(- 


l + 2w 2 )(9 + 4w 2 )(l-5w 2 )(- 
758804^^6 


l + 2co 2 )(9 + 4a> 2 ) 




5- 2 )(- 


-1 + 2cj 2 )(9 + 4w 2 )(1-5w 2 )(- 
13040 


-l + 2co 2 )(9 + 4u; 2 ) 


1 2(1 


-5o, 2 )( 


-1 + 2w 2 )(9 + 4cj 2 )(1-5cj 2 )( 


-l + 2co 2 )(9 + 4w 2 ) 
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1 j . f 58477 89216 7872 



294912 1 1(1-2^2)2 (9-14^-8^) (-l + 2u; 2 ) 2 
33456 



+ 



(9 - 14c 2 - 8w|) 



} 



. o o f 5864788 186165 1885814784 

+ 2^1^ n 



\(l-2w 2 ) (-l + 2c 2 ) 2 (4 w 2 -c 2 2 ) + (c 2 -4c 2 )(9-14c 2 -8u4) 
18210816 1 111689728^| 1 



+- 



(1 - 7c 2 + 10wf)(l - 7w| + 10w|) J (1 - 7w 2 + 10wJ)(l - 7w 2 + 10w: 



D5 = ' J-„ 24 f,„ + 6426 



49152 | \ (l-2w 2 ) 2 (-9 + 14c 2 + 8uj ) 

30450688 

~~ (-1 + 5c 2 ) 2 (9 - 14c 2 - 8cf)(-l + 5c 2 ) 2 (9 - 14w 2 - 8w|) 
2 2 f 90048 139298 39249 



(l-2c 2 ) 2 (l-2c 2 ) 2 (l-2c 2 ) 2 (4c 2 -c 2 ) 



447897600 952565760 

+ 7 - . » o7^ . ; ^ + 



(-l + 2c 2 (9 + 4c 2 ) (9 - 14c 2 - 8c|) 
6276089856 



(1- 


5c 2 ) 2 (-9 + 14c 2 + 8cf)(l - 5c 2 ) 2 (- 


-9 + 14c 2 + 8c|) 




594542592 ) 






-4c 2 ) (-9 + 14^ + 8^) J 






3159788544wfw| 




+ U- 


5c 2 ) 2 (-9 + 14c 2 + 8cf)(l - 5c 2 ) 2 (- 


-9 + 14c 2 + 8cf) 




49312045056^ 




+ d- 


5c 2 ) 2 (-9 + 14c 2 + 8cf)(l - 5c 2 2 ) 2 (- 


-9 + 14c 2 + 8cf) 




3734949888c?cf 





(1 - 5c 2 ) 2 (-9 + 14c 2 + 8cf)(l - 5c 2 ) 2 (-9 + 14c 2 + 8c|) 



^ 6 82944V3{ 29889 { (-1 + 2c 2 ) 2 + (9 - 14c 2 - 8c 2 4 ) } 

. 2 2 f 738 93899 91445760 



(l-2c 2 ) 2 (-l + 2c 2 ) 2 (c 2 -4c 2 )(9-14c 2 -8ct) 
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D 7 = ,^,127/™*+. ^ 



110592^1 I (-l + 2w?) 2 (w 2 -4u; 2 )(-l + 2 w 2 ) 2 
758086 5904 1 



(w 2 -4w 2 )(-l + 2^ 2 ) (9-14w|-8w|)J 

_ 2 of 492 370964 58653 

+2^| J- + 



+ 



l-2w 2 ) 2 (-l + 2w 2 ) (4w 2 -w 2 )(-l + 2u; 2 ) 2 
13893120 116702208 



(9 + 4cj 2 ) 2 (1 - 2cj|) (co 2 - 4w 2 )(-l + 2co 2 ) 2 (9 + 4w 2 ) 2 



103680 870912 



(9 + 4w 2 )(l-2^ 2 ) 2 (a; 2 -4w 2 )(-l + 2a ) 2 ) 2 (9 + 4u; 2 ) 2 
62519040 246177792 



(-9 + 14w 2 + 8a4) (w 2 -4w 2 )(9 + 4co 2 ) 2 
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